Crossover from adiabatic to sudden interaction quench in a Luttinger liquid 
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Motivated by recent experiments on interacting cold atoms, we analyze interaction quenches in 
Luttinger liquids (LL), where the interaction is ramped from zero to a finite value within a finite 
time. The fermionic single particle density matrix reveals several regions of spatial and temporal 
coordinates relative to the quench time, termed as Fermi liquid, sudden quench LL, adiabatic LL 
regimes, and a LL regime with time dependent exponent. The various regimes can also be observed 
in the momentum distribution of the fermions, directly accessible through time of fiight experiments. 
Most of our results apply to arbitrary quench protocols. 
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Non-equilibrium dynamics and strong-correlation phe- 
nomena in quantum many body systems are topics at the 
forefront of contemporary physics. When these two fields 
are combined, namely when strongly correlated systems 
are driven out-of-equilibrium, we face a real challenge. 
Experimental advances on ultracold atoms [l[ have made 
the time dependent evolution and detection of quantum 
many-body systems possible, and in particular, quantum 
quenching the interactions by means of a Feshbach reso- 
nance or time dependent lattice parameters has triggered 
enormous theoretical j^-Q and experimental activ- 
ity. 

Luttinger liquids (LL) are ubiquitous as effective low- 
energy descriptions of gapless phases in various one- 
dimensional (ID) interacting systems [13, Hi- In ID 
fermions, e.g.. Landau's Fermi liquid description breaks 
down for any finite interaction, and the low-energy 
physics is described by bosonic collective modes with lin- 
ear dispersion, and is characterized by anomalous non- 
integer power-law dependences of correlation functions. 
The LL similarly arises as the low-energy description of 
interacting ID bosons or that of spin chains [loj . 

A LL could also be driven out of equilibrium through 
transport 12, [l^, e.g., but here we shall concentrate 
on time dependent changes of the interaction parameter, 
which is of particular relevance for cold atomic systems. 
Sudden quenches (SQ) of the interaction in LLs have 
been considered recently by several authors 14-l3|, and 
the idea of the other extreme limit of adiabatic param- 
eter ramps is often invoked, too. However, experimen- 
tal ramps cannot take infinite time, and are not instan- 
taneous, either. Here we study, how a nonzero quench 
time, T 7^ 0, influences the final state of the system af- 
ter a quantum quench. As we show, a finite r leads to 
'heating' effects, and generates excitations in the final 
state. Moreover, it amounts in the appearance of addi- 
tional energy (^ I/t) and corresponding length scales: 
while in certain space-time regions the system reveals 
universal near-equilibrium (adiabatic) correlations [H , in 



other regimes renormalized Fermi liquid (FL) or sudden 
quench (SQ) behavior is found. 

Thus motivated, let us study the LL Hamiltonian[l(^ 

H = ^a.(g)6t6^ + ^[6,6_, + 5+5t J, (1) 

95^0 

with Lo{q) = v\q\ [v being the bare "sound velocity"), 
and 6j the creation operator of a bosonic density wave. 
The interaction g is changed from zero to a nonzero 
value within a quench time r, g{q, t) = g2((l)\q\Q{t), with 
Q{t) encoding the explicit quench protocol, and satisfy- 
ing Q{t > r) = 1 and Q(t < 0) = [H. In particular, 
for a linear quench, Q{t) = tQ{t{T — t)) /t + Q{t — t) with 
Q{t) the Heaviside functions. This setting is very gen- 
eral, and equally applies to switching on the interaction 
in a spinless fermion system, quenching away from the 
Tonks- Girardeau limit of a ID Bose gas Q, or turning 
on the Z-component of the interaction in an XXZ spin 
chain [loj. To make contact with previous work [l^ - [l6| . 
here we focus on fermionic correlators in detail, but most 
of our results can be trivially generalized to interacting 
bosonic systems [18 1. 

We describe time-evolution using the Heisenberg equa- 
tion of motion, leading to 



idtbq = [bq, H] = uj{q) bq + g{q, t) btq 



(2) 



and similarly, idtb~^„ = —uj{q)b't„ — g{q,t)bq. Solutions 



of these are of the form 



bq{t)=u{q,t) bq{0)+V*{q,t) btqiO) 



(3) 



where all the time dependence is carried by the prefac- 
tors, u{q,t) and v{q,t), and the operators on the r.h.s. 
refer to non-interacting bosons before the quench. All 
expectation values are thus taken in terms of the ini- 
tial density matrix of the latter (or vacuum at T = 
0). The bosonic nature of the quasiparticles requires 



2 



- = 1. From Eqs. ©-(O, wc obtain 



idt 



u{q,t) 




_ v{q,t) 





-g{q,t) -^{q) 



u{q,t) 
v{q,t) 



(4) 



with the initial condition u{q,0) = 1, v{q,0) = 0. Since 
both Ltj{q) and g{q,t) are even functions of q, u{q,t) and 
v{q,t) must be so, too. By Eq. all time depen- 
dence has been transferred to the Bogoliubov coefficients, 
and therefore expectation values of the time dependent 
bosonic modes and non-equilibrium dynamics are calcu- 
lable using standard techniques developed for equilibrium 
[lo| , once the solutions of Eq. ([4]) are known. 

Before discussing a continuous quench, let's see how 
limiting cases are recovered from Eq. Q . The adiabatic 
limit follows from replacing g{q, t) with its time indepen- 
dent final value, and looking for the stationary solutions 
of Eq. (|4]) at a given energy while ignoring the initial 
conditions. The SQ limit requires only the replacement 
of g{q, t) by its final, time independent value, and solving 
the resulting linear differential equation with the initial 
conditions satisfied. 

With a linear quench, Eq. ffl realizes the non- 
hermitian Landau-Zener model [19[, which can be solved 
exactly in terms of the parabolic cylinder function. How- 
ever, the exact solution does not yield an immediate and 
transparent physical picture. Therefore, to obtain more 
insight, we assume that g{q,t) ^ w(q) (i.e., g2{q) ^ v) 
for all t and q, and solve Eq. (|4]) perturbatively in 
the interaction. To lowest order in g2{Q), we obtain 
u{q,t) « e'}q>{—iuj{q)t) and 



v{q,t>0) 



II dt'g{q,t')exp{iuj{q){t-2t')) . (5) 
/o 



Higher order corrections to u{q, t) and v{q, t) are of order 
and 52(9)1 respectively. We have also checked nu- 
merically that Eq. ([5]) is indeed applicable for any t and 
T, as long as g2{q) ^ v |3J|. In the SQ (t 0) and 



adiabatic (r 00) limits we obtain 



.(<Z,t>r)«.^x 
Zv 



2isv[\{u){q)t) for r ^ 0, 
- exp(— ia;(g)i) for r — > 00, 



reproducing to lowest order in 52(9) the SQ results 14, 
16 1 and the equilibrium Bogoliubov transformation 10, 
2l| . respectively. 

We are now in position to obtain information about 
physical observables. We start with the evolution of the 
total energy of the system. We take the energy of our 
initial vacuum state to be zero. In the fermionic setting, 
this corresponds to measuring the energy with respect to 
the energy of the non-interacting Fermi sea. The expec- 
tation value of Eq. ([T]) is then evaluated in the Heisen- 
berg picture, where the expectation value is taken with 
the non-interacting ground state, and v{q,t) and u{q,t) 



as obtained from Eq. ^ keeping track of the time evo- 
lution of the system. Wc thus obtain for {H{t)) after the 
quench (t > 0) 

(H) =J2^iqhB{q) + {2nB{q) + l)Mv*{q,t)dtv{q,t)] , 

with nsiq) — l/(cxp(aj((7)/r) — 1) the Bosc function. 
The expression above is time independent for i > r, as 
expected. At T = and t > t, and an interaction of 
finite range, (72(9) — ff2 cxp(— _Rokl/2), we obtain 



(H) = En 




(7) 



for a linear quench. Here we introduced the microscopic 
time scale, tq = Ro/2v, and Egg ~ — Lg|/47rui?Q is the 
adiabatic ground state energy shift to lowest order in g2, 
with L the system size. The second term corresponds to 
quasiparticle excitations resulting from the finite quench 
speed. In the SQ limit, r ^ tq, the energy of the system 
is only slightly shiftedjli], (H) = Egs{T/Tof/2. This 
holds true for a general quench, i.e. (iJ) ~ (t/tq)"^ when 
T — > with a quench dependent coefficient. In the adia- 
batic limit, T ^ To, on the other hand, the excess en- 
ergy (or "heating") vanishes as —2Egs\n{T/To) Tq/t^ 
in accord with the so-called analytic response of Ref. 
[2^. This remains valid for general smooth quenches 
displaying kink(s) (discontinuity in the derivative) and 
bounded dtQ(t). Smooth quenches without kinks but 
with bounded dtQ{t) produce also a universal decay as 
^ 1/t^, while the r dependence of the heating becomes 
non- universal for protocols with a diverging dtQ(t)\24\. 
The crossover between the SQ and adiabatic limits oc- 
curs when T ^ Tq, which typically translates to r 1/J 
in an optical lattice, with J the hopping integral in the 
underlying microscopic lattice Hamiltonian. 

In the fermionic context, the structure of the non- 
equilibrium dynamics can be well demonstrated by means 
of the fermionic one-particle density matrix. Since the 
fermion field decomposes to right-going and a left-going 



(•g^ parts, 'I'(a;) 



^r{x) + e~''=-^^*i(x), it is enough to 



concentrate on the right-going part of the density matrix, 
Gr{x,t) = (*+(a;,t)*,(0,i)) , 

describing excitations around the right Fermi momen- 
tum, k sa kp. The right-going field, ^'r(a;), can be ex- 
pressed in terms of the LL bosons as [loj 



Vr 



: exp (i4>r{x)) 



(8) 



where i]r denotes the Klein factor, and (prix) = 
J2q>n \/WM^e'9^-"l''l/2&g _^ Following standard 
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steps U I , we obtain then 



g>0 ' 



Grix,t) = G°{x)cxp 

X [nB{q) + Hq,tf\{2nB{q) + l)]), (9) 



where G°(a;) = i/(27r(.T + ia)) denotes the free fermion 
propagator, with a an ultraviolet regulator. At T = 0, 
the Bose functions vanish, and only v{q,t), i.e. the mix- 
ing between bq and determines the dynamics. 

Let us first discuss the properties of Gr{x,t) long af- 
ter the quench, t ^ t. In this limit, we can show (25j 
that, independently of the quench protocol, Q{t), the 
onc-particle density matrix exhibits universal properties. 



Grix,t) 



A{r/To) 



Rn 



min{|a;|, 2vt} 



Rq 



7SQ 

for |x| > 2vT, 

for < 2ut, 
(10) 



where 7sq = g^/v^ + ... and 7ad = g|/2u^ + . . . de- 
note the perturbative sudden quench and adiabatic ex- 
ponents, respectively. The prefactor A{t/tq) depends 
on the speed of the quench: For a sudden quench it is 
A{t <C To) ^ 1, while for slower quenches A{t > tq) ^ 

Thus even for t — > 00, instead of one single power- 
law, Gr interpolates between the SQ and adiabatic limits. 
This is shown in Fig.[T]for a linear quench j25|. Physically, 
it is easy to understand the cross-over behavior observed 
in Gr : a finite-time quench is experienced by slow excita- 
tions of energy Lj{q) < 1/r as a sudden change, while fast 
excitations with Lj{q) > 1/r can adjust to the change in 
the interaction strength adiabatically. Since high (small) 
energy excitations determine the short (long) distance 
correlations, the tail of Gr is governed by the SQ expo- 
nent [l3|, while the short distance behavior is described 
by the adiabatic exponent. It is remarkable that for slow 
enough quenches, r 3> tq, the quench time manifests it- 
self explicitly through an adiabatically enhanced prefac- 
tor A ^ (t/tq)'^'^'^ of the asymptotic tail as also shown in 
Fig. [TJ Thus while the spatial decay of Eq. (fTOl) contains 
the SQ exponent, its r dependence reveals the adiabatic 
LL exponent. For a finite t ^ t but 2vt ^ |a;|, Gr{x,t) 
decays asymptotically as iZ(t)/2TTX, with a finite quasi- 
particle weight. 



Z(t»r,ro)^A(r/ro) (^) 



7SQ 



(11) 



Thus the exponent observed in Z{t) is identically 7sq for 
t 2> T, but the finite quench time amounts in a quasi- 
particle weight increased by a factor, A ^ (t/tq)^'"^ for 
T ^ tq. Although these results were obtained perturba- 
tively, they carry over to the non-pcrturbativc limit, too. 
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FIG. 1: (Color online) Left: the long time {t <^ t ^ 00), 
steady state limit of the one-particle density matrix is plotted 
on loglog scale for a linear quench and 2(j2 = w as a function 
of |a;|, exhibiting the crossover from adiabatic behavior (lower 
line in Eq. (Iiop l at small to SQ behavior with adiabatic 
enhancement (upper line in Eq. (|10|) ') at large |a-|. The curves 
are plotted for r/2ro = 10*, 4 x 10'* and 16 x 10"' from bottom 
to top in both panels. Right: Landau's quasiparticle weight, 
Z{t) is plotted on a loglog scale as a function of t, bridging 
between the weakly interacting Fermi liquid to strongly sup- 
pressed Z <^ 1 with adiabatic enhancement (Eq. (|lip ). Inset: 
the linear quench protocol is shown. 



with the only exception that the exponents 7sq and 7ad 
must be replaced by their exact value. 

All these spatial features appear also in the time- 
dependent momentum distribution of the fermions, 
n(k,t), directly measurable through time of flight exper- 
iments. In particular, at T = and finite t, n(k,t) 
exhibits a jump of size ~ Z{t) at k = kp, while it ap- 
proximately scales for |fc| 3> l/2vt as 



n{k) 



-sign(fc) X 



A(r/To) |fci?or^«, 1^1 « 

lfci?^r■'^ 1^1 » 



1 

2VT ' 

1 

2vT ' 

(12) 



ioT k = k — kp, \k\ -^i kp, and t ^ t. Thus the time scale 
of the quench is also imprinted in the momentum distri- 
bution, which also shows a cross-over behavior between 
the SQ and the adiabatic limits. For adiabatic quenches, 
T — > 00, we recover the equilibrium LL exponent, while 
close to kp, the momentum distribution is enhanced by 
a factor A{t/tq) compared to the SQ behavior l4 llq. 

The above analysis can be extended to the short time 
region, t ^ r, where the behavior found depends explic- 
itly on the quench protocol 25j as 



Gr{x,t) - G°(x) 



Rn 



min{|a;|, 2vt} 



lit) 



(13) 



where ^[t) = g|(5'^(i)/2w^ -I- For short distances, 

|a;| ^ 2vt, the spatial correlations decay with a time- 
dependent exponent, and this region can thus be charac- 
terized as a weakly interacting LL (t-LL). For \x\ ^ 2vt, 
on the other hand, similar to t 3> t, correlations re- 
main almost unaffected by interaction, and a Fermi liq- 
uid regime is found. For i ^ tq, .^(t) ~ 1 as in the initial 
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Fermi gas, but for t ^ tq we recover a Fermi liquid be- 
havior with a reduced quasiparticle weight as 



feature tunable interactions, essential to address quench 
dynamics. Amon g th ese, ID configurations have been 



Z{to « < « r) ~ 



realized using '*°K[26|, l27|, ''Li [291, and the momentum 



(14) 




2VT 



FIG. 2: (Color online) The schematic universal spatial- 
temporal characteristics of a quenched LL, with the bound- 
aries denoting crossovers. In the adiabatic LL regime, the LL 
exponent of the final state, 7ad governs spatial correlations, 
while in the SQ LL region, correlations decay with the SQ ex- 
ponent, 7sQ > 7ad. Correlations are adiabatically increased 
by an amplitude, A ~ (r/ro)"'-'"' in the shaded region. In the 
Fermi liquid region a time dependent quasiparticle residue is 
found (see Eq. (|14|) ). while in the time-dependent LL (t-LL) 
region a quench protocol-dependent weakly interacting LL is 
found with a time dependent exponent. Eg. (|13|l . The dashed 
line denotes = 2vt, i.e. the light-cone For r <C tq, the 
SQ physics of Ref. 

dominates everywhere. 

The quasiparticle weight thus slowly decreases during 
the quench, and excitations remain similar to those in 
the initial Fermi gas with a reduced weight {Z < 1) for 
t < T. After the quench, t > t, the quasiparticle weight 
continues to decrease as a power-law, and it resembles 
to an interacting (heavy) Fermi liquid with 1/x spatial 
decay and Z <^ 1 quasiparticle residue, Eq. pT|) . as was 
also found for sudden quenches [H, [3l ■ This situation 
is shown in detail in Fig. [TJ where Z{t) is plotted for 
the special case of a linear quench jisl]. Our results are 
summarized in Fig. [2] 

The main effect of finite temperatures (T > 0) in Eq. 
^ is to introduce a new time or length scale, 1/T or 
v/T, respectively, above which the power law behaviour 
of LL disappears, and gives way to exponentially sup- 
pressed behaviour as exp(— crmax[|a;|/w, t]) with c some 
constant. Our findings in Fig. [5] survive in the region 
it,T,\x\/v) < 1/T. 

Ultracold fermionic gases have been realized using sev- 
eral atoms such as ''OKjii^, ^Lijli, i"Yb-"3Yb[3|, 
i^^Dyim and ^'^Srjiil, and temperatures well in the 
quantum degeneracy regime were reached (T < 0.1 Ep, 
with Ep the Fermi energy). All these atomic systems 



distribution has been measured in time of flight (ToF) 
experiments in 2D[331 and 3D[2^,[3^ Fermi gases. There- 
fore, by applying ToF imaging or momentum resolved rf 
spectroscopy [26|, the observation of the momentum dis- 
tribution of Eq. (fT2|) is within reach for ID fermions [2^. 
Furthermore, the specific momentum distribution of a LL 
has already been observed in the Tonks- Girardeau limit 
of ID Bose systems 0, which exhibit fermionic proper- 
ties in this strongly interacting regime, and reveal after 
an interaction quench features similar to the ones found 
for fermions [18|. 

In summary, we have studied continuous interaction 
quenches in LL, bridging smoothly between the SQ and 
adiabatic limits. The resulting dynamics is largely infiu- 
enced by the finite quench time for fermions, and in par- 
ticular, the momentum distribution exhibits a crossover 
from the adiabatic LL to that of the SQ with the extra 
adiabatic enhancement factor (t/tq )''"''', revealing both 
the equilibrium and SQ LL exponents. A finite quasipar- 
ticle residue is retained during the quench, refiecting the 
Fermi gas nature of the initial state, getting suppressed 
gradually after the quench. The variety of quench in- 
duced phases offers a unique opportunity to design low 
dimensional correlated states on demand. 
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SUPPLEMENTARY ONLINE MATERIAL FOR "CROSSOVER FROM ADIABATIC TO SUDDEN 
INTERACTION QUENCH IN A LUTTINGER LIQUID" 

EXACT EVALUATION OF THE BOSONIC CORRELATOR OF EQ. (9) IN THE MAIN TEXT FOR A 

LINEAR QUENCH 



For a linear quench, we can integrate Eq. (5) in the main text to obtain 

v{q,t) = 1^^^ [sin{u{q)t) - uj{q)tcx^^{-iu{q)t)] 



for < i < T, and 



2Lj2(g)r 



viq,t) = ^f^l^*^!^^ [exp(za;(g)(t - 2r)) - exp(za;(g)f) + 2ta;(g)T exp(~za;(g)t)] 



(15) 



(16) 



for t > T. 

The exponent in the one-particle density matrix in Eq. (9) in the main text is then evaluated in closed form for 
< > T, T = and L cx) as 



q>0 



(f) = -§ihir,x,Ro) + l2it,T,x,Ro)), 



where 



(17) 



Ii{t,x,Ro) ^ dq 



exp(-i?o(7) . 2 f<l^\ 2 



q3y2^2 



(^) [sin\ujiq)r) + (u:iq)rr]. 



l2it,T,x,Ro) = dq 



exp(-i?oq) . 2 



q^vT 



sm 



(^) [- sin(2cj(g)t) + sin(2t^(g)(i - r))] , 



(18) 
(19) 



which are evaluated as 

Ii{t,x,Rq) 



32(ut)2 



{ln(2irt;r + i?o + isx){sx — iRo + r2vT)^} 



2 {hi(i?o + isx)[{Ro + isxf + {^vrf] + \n{2isvT + i?o)(i?o + 'iisvTf] - 41n(i?o)[-Ro + (2wt)^ 



s=±l 



(20) 



and 



hit, T, x, Ro) = i In {Rl + 4i>2(t - tf) - ^ ^li^ (i?o + (2wt - 2vt - sx)^ 

s=±l 



E 

s=±l 



2{2ivt - sRo) In 



sRq + 2iv{T - t) 
sRq — 2ivt 



(rRo — 2ivt + isx) In 



r=±l 



rRo + i{2vT - 2vt + sx)^ 
ri?o ~ i{2vt — sx) 



(21) 
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The < t < T correlator is obtained as 

-E (jz) (f ) l''^'^'^)'' = " iIiit^^^Ro)+Ht,t,x,Ro)) . (22) 

Expanding these in various hmits are used to obtain the resuhs cited in the paper, and their general {t,T, x, Ro) 
dependence is used to generate Fig. 1. in the main text. 



EVALUATION OF THE ASYMPTOTICS OF BOSONIC CORRELATOR OF EQ. (9) IN THE MAIN 

TEXT FOR GENERAL QUENCH PROTOCOL 

In this section, we demonstrate that the asymptotic behavior what we obtained for a linear quench is universal, 
and arbitrary quench protocols lead to the same behaviour. For arbitrary quench protocol Q{t), the exponent in the 
one-particle fcrmionic density matrix of Eq. (9) in the main text at T = and L ^ oo can be rewritten as 



/ = - 



i:(|)-'"'(f)w-')'i-3j/*/* 

g>0 ^ n n 



Q{ti)Qit2)dtAJiti~h), 



where 



fit) = In 1 + 



(Ro - 2ivty 



After partial integrations, it takes the form 



(23) 



(24) 



I = - 



4^2 



Q\t)f{0)-2Q{t) / dhQ'{ti)Rcf{ti-t)+ / dtiQ'iti) / rfi2Q'(i2)./(ti - i2) 



(25) 



Here, by using Q{t > r) = 1, the upper limit of integration would reduce to min{t, r}. However, our considerations 
remain valid for smooth quench functions as well, reaching 1 only asymptotically, i.e Q{t ^ t) — > 1. Let us first 
consider the properties of the steady state in the limit of t (r, x/u), when the middle term in Eq. does not 
contribute, and Q{t ^ r) = 1. For t ^ \x/v\ ^ r. the first term yields 2 ln(|a;|/i?o)j while the last integral produces 
similar spatial decay and the adiabatic enhancement as 2 ln(|a;|/2wr). The exponent is 



for t > \x/v\ > T. 



VR^t 

In the t ^ T ^ limit, the last term also vanishes, and we are left with the adiabatic exponent 

X 



I 



2^2 



In 



Ro 



for t > T > \x/v\. 



(26) 



(27) 



The other limit of interest is \x\ ^ {vt, vt), when the exponent simplifies in Eq. with /(t) replaced by /i(t) = 
— 21n(i?o ^ 2ivt). The first term always only contributes with a constant, — 21n(i?o)- In the limit of |x| ^ vt ^ vt, 
Q{t) = 1, the third term becomes independent of both x and t and gives rise to the adiabatic enhancement factor as 
— 21n(T), and only the second term determines the temporal decay. The exponent is obtained as 



9l 



In 



/ToT 



for \x\ ^ vt VT 



(28) 



In the limit of |a:;| 3> vt 3> vt, i.e. during the quench, the time dependence of Q(t) is essential. The second term 
yields AQ^{t) ln(t) to leading order, while the third term gives —1/2 times the second term. Altogether, the exponent 
reads as 



/ = -^Q^{t) In ( — ) for > wr > vt 



2w2 



To 



(29) 



These agree with the asymptotic expansion of the exact res ults for a linear quench. 



